Semiclassical theory of shot noise in disordered SN contacts 
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We present a semiclassical theory of shot noise in diffusive superconductor - normal metal contacts. 
At subgap voltages, we reproduce the doubling of shot noise with respect to conventional normal- 
metal contacts, which is interpreted in terms of an energy balance of electrons. Above the gap, 
the voltage dependence of the noise crosses over to the standard one with a voltage-independent 
excess noise. The semiclassical description of noise leads to correlations between currents at different 
electrodes of multiterminal SN contacts which are always of fermionic type, i.e. negative. Using 
a quantum extension of the Boltzmann - Langevin method, we reproduce the peculiarity of noise 
at the Josephson frequency and obtain an analytical frequency dependence of noise at above-gap 
voltages. 



PACS numbers: 72.70.+m, 74.40+k, 74.50+r 

In recent years, the noise properties of hybrid con- 
tacts involving superconducting (S) and normal (N) met- 
als attracted considerable attention.^ Theoretical work 
in this direction was pioneered by Khlusjj who found 
that Andreev reflectionjj in which electrons incoming 
from normal metals are reflected from the SN interface 
as holes, play a key role in the shot noise of SN struc- 
tures. For clean SN contacts, he predicted that the shot 
noise vanishes ^at subgap voltages eV < A. de Jong 
and Beenakkeia addressed dirty SN contacts and found 
that the shot noise at subgap voltages is doubled with 
respect to its value in a normal contact with the same 
resistance. This doubling mL shot noise has already been 
experimentally confirmed. BU The distribution function of 
current fluctuations was calculated and it was found that 
at subgap voltages, it describes independent transfers 
discrete charge 2e through the contact .0 Lesovik et al 
found that the frequency dependence of the shot noise ex- 
hibits a peculiarity at the Josephson frequency to — 2eV 
instead of u> = eV in normal contacts. More recently, the 
same authors obtained additional singularities in the fr< 
quency dependence of the noise at above-gap voltages 
These results could lead to the impression that the mo- 
tion of electrons and holes is correlated not only at the 
SN interface, but that the current is carried also through 
the normal part of the contact in portions of 2e instead 
of e. Recent findings that the correlations of currents at 
two normal, con tacts attached to a superconductor may 
be positive JlJ'tll as in the case of Bose statistics seemed 
to reinforce such a view. On the other hand, from these 
works it is not clear whether the above results require 
phase coherence of electrons nor if there are restrictions 
for the contact length. 

Recently, Gramespacher and one of the authorslll dis- 
cussed the fluctuations in multiterminal SN contacts 
using clcctronlike and holelike distribution functions, 
which were expressed in terms of quantum-mechanical 



injectivitiesEj of the terminals. It was shown there, in 
particular, that positive correlations vanish in the limit of 
a large channel number. In this paper, we present a semi- 
classical theory of shot noise in disordered SN contacts. 
We show that the doubling of shot noise may be obtained 
within a simple Boltzmann - Langevin approacfo pro- 
vided that the appropriate boundary conditions for the 
average distribution function are used. We relate the in- 
crease in the shot noise to the peculiar energy transport 
through an SN interface and show that this effect is stable 
with respect to phase-breaking. From the semiclassical 
nature of electron transport in the N region, it imme- 
diately follows that the correlations between currents in 
different electrodes of a multiterminal diffusive SN con- 
tact are always negative. We are also able to explain 
the peculiarities of the frequency dependence of noise in 
terms of the shape of the electron distribution function. 



Consider a narrow normal-metal microbridge connect- 
ing massive normal and superconducting electrodes. The 
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FIG. 1. The non-equilibrium electron distribution func- 
tion in the normal part of an NS contact at voltage 
eV = I.4A. 

elastic mean free path I of electrons in the microbridge 
is assumed to be short and the length of the contact L is 
assumed to be much larger than (HD/A) 1 / 2 , where A is 
the superconducting gap and D is the diffusion coefficient 
in the normal metal. The applied voltage or temperature 
are considered to be much larger than the Thouless en- 
ergy e t — hD/L 2 . Under these conditions, it is possible 
to neglect the penetration of the condensate into the mi- 
crobridge and consider it just as a normal metal . with a 
nonequilibrium distribution function of electrons .r 5 H 16 l 

Let the x axis be directed along the contact, x — cor- 
responding to the superconducting electrode, and x = L, 
to the normal-metal one. The electron energy s in the 
normal metal is measured from the Fermi level of the 
superconductor. Introduce an energy-dependent coeffi- 
cient W of normal transmission of an electron from the 
contact through the NS interface without being Andreev 
reflected, so that W{e) = if |e| < A and W(e) -> 1, for 
|e| » A. 

Since the electrons in the superconducting electrode 
are described by an equilibrium Fermi distribution func- 
tion /o(e), the distribution function of electrons moving 
into the cawtact from the SN interface may be written in 
the fornJliy 

f(e, v x ) = Wfo(e) + (1-W)[1- f(-s, -v x )}, (1) 

for e > and v x > 0, and the distribution function of 
holes moving into the contact from the interface may be 
written in the form 

1 - /(e, v x ) = W[l - fa(e)] + (1 - W)f(-e, -v x ), (2) 

for e < and v x > 0. Since in the bulk of the contact 
the electron motion is diffusive, we must now obtain a 
suitable boundary condition at the NS interface for the 
diffusion equation. 

Introduce the symmetric and antisymmetric parts of / 
with respect to v x : 

/.(e) = (1/2) [f(e,v x )+f(£,-v x )}, 

/„(£) = (1/2) [/(£,<) -/(£,-<)]. (3) 



[ + Me-eV) + \(l-%)Me + eV), 
{ (l-f)/o(e) + f/o(e-en 



At subgap voltages, the shape of the distribution func- 
tion is exactly the same as if the contact were extended 
in negative direction frorruO to — L and the voltage — V 
were applied at this point .□ The distribution function for 



Inside the contact, the distribution function obeys the 
diffusion equation and may be represented as a sum of 
an isotropic part /(e, x) and a small anisotropic part 

fi = -(v x /v F )(l/L)[f(£,L)-f(£,0)]. (4) 

Following Kupriyanov and Lukichev,ll§ we equate now 
/(e,0) to f„{e) and fx to (v x /v F )f a (£), Introducing the 
distribution function of electrons in the normal electrode 
f n (s), one may express f a in terms of f s and /„: 

f a (£) = (l/L)[f s ( S ) -/«(£)]. (5) 

For the sake of dcfiniteness, assume now e > 0. By 
substituting Eqs. (S) into Eqs. (0) for f(£,v x ) and (g) 
for 1 — /(— e, v x ) and excluding f a from them by means 
of (||), one obtains a closed set of equations for f s (e) and 
fs(-e): 

(1 + l/L) f s (e) + (l- W) (1 - l/L) f s (-£) = 1 - Wfo(-e) 

+ (l/L)f n (£) - (l/L)(l - W)f n {-£) (6) 

and the equation obtained from it by interchanging e and 
— e. The solution of the set (||) assumes different forms 
depending on the relationship between W and l/L. In 
the subgap region |e| < A, W — and we obtain 

fs{e) = \[l + Ue)-fn{-e)] (7) 

for both signs of e. Suppose that the normal electrode 
has a potential V. Then / n (e) = fo{e — eV), and inside 
the gap, 

fs{e)= 1 -[h{e -eV) + f Q (£ + eV)], \e\ < A. (8) 

Outside the subgap region, W ^ 0, and l/L may be ne- 
glected with respect to it virtually at all energies. In this 
case, one easily obtains 

f.(e) = f (s), \e\ > A. (9) 

Inside the contact, the distribution function is obtained 
by solving the diffusion equation with the boundary con- 
ditions and @. As a result, one obtains: 

\e\<A 

(10) 

\£\>A. 



an above-gap voltage is shown in Fig. 1. Note that at 
finite voltages, /(e) is discontinuous at |e| = A even at 
nonzero temperature. 

Consider now the noise of the contact. Since all the 
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voltage drop takes place inside the normal-metal micro- 
bridge and there is no voltage drop across the SN inter- 
face, it is reasonable to assume that all the noise is pro- 
duced by random impurity scattering in the microbridge 
only. In this case, we can use the welLknown Langevin 
equations for the current fluctuations £3 

5i = -D^-5p-a^-84> + 5r\ (11) 
or or 

where D is the diffusion coefficient, a is the electric con- 
ductivity, Sp(r) is the charge-density fluctuation, 6<j>(r) 
is the local fluctuation of the electric potential, and the 
correlator of extraneous currents 8'f xt is given by 

(«Sf*(n)«r*fa)><- = 4<r<W£(ri - r 2 )Ttf(n), 



T N (r) = J de/(e,r)[l-/(e,r)]. (12) 

At sufficiently low frequencies, the fluctuation of the to- 
tal current <5j may be considered as constant along the 
contact, and Eqn. flll] ) may be integrated over its length. 
As a result, the gradient terms in ( pj] ) drop out because 
of the boundary conditions Sp — and Scf) = 0, and one 
arrives at the standard, expression for the noise in disor- 
dered metal contacts:EH 

^ = |/t Tjv(2;) ' (13) 

where R is the normal resistance of the contact. Substi- 
tuting now the distribution function ([To]) into (|T3|) and 
performing the integration, one obtains that 
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exp(A/T) + exp(eV/T) 



exp(A/T) + cxp(- e y/T) 



(14) 



At zero voltage, this expression reduces to the Nyquist 
formula Sj — AT/ R. At zero temperature yet finite volt- 
age, Eqn. (|I4|) takes a very simple form: 



o 4e|V| 
Sl= 3~R- 



~9(e\V\ 



A)- 



e\V\ 



R 



(15) 



The shot noise in diffusive SN contacts is doubled at volt- 
ages e\V\ < A as compared to the same contacts between 
normal metals. The reasons for the shot-noise doubling 
can be understood in terms of the noise temperature of 
electrons Xjy and the energy transport. The electrons 
in the contact acquire additional energy because of their 
acceleration by the field. If both electrodes are made 
of normal metal, the excess energy flows into both elec- 
trodes thus effectively cooling the contact and decreas- 
ing Tjy. However if one of the electrodes is supercon- 
ducting, the SN. boundary forbids heat transfer in the 
subgap region. E3 The absence of heat transfer through 
an SN boundary has been known for many years, and 
this is precisely what motivated the pioneering paper by 
AndreevE Because there is only one energy drain now 
instead of two, naturally T/v is higher than in the case of 
a normal contact. 

The above semiclassical derivation shows that the dou- 
bling of shot noise is insensitive to phase-breaking. In the 
case of strong electron-electron scattering, the effective- 
temperature approximation may be used at eV <C A 
and precisely the same heat-balance equationJjoij the ef- 
fective temperature T e as in normal contactsEjtfl can be 
written except that the heat-absorption boundary con- 
dition T e = T at the SN interface should be replaced 
by the zero heat-flux condition dT e /dx = 0. It is eas- 
ily seen that this should double the shot noise making it 



2 x (v / 3/2)e/ = \fZeI and thus increasing it above the 
nonintci ja tj ti pg-clcctron value precisely as for a normal 
contact.EjEa Much like in the case of normal contacts, the 
noise should be suppressed by energy relaxation caused 
e. g. by phonon emission. 

Using the semiclassical approach, Sukhorukov and 
Lossc3 recently investigated the correlations between 
currents in different contacts of multiterminal diffusive 
normal-metal systems. Here we extend theic discussion 



to the case of multiterminal SN systems.cSl It-, is con- 
venient to introduce characteristic potentialsE3 </> n (r), 
n = I, . . . , N associated whith each contact n which equal 
1 at contact n and are zero at all other electrodes obeying 
the equation 



V(aV0„) = 



(16) 



in the bulk of the contact. In terms of these potentials, 
the normal-state conductance matrix of the contact may 
be written in the form 



G m n = ~ I dS m ■ <7<fi m V(f) n , 



(17) 



where S m is the interface of the system with electrode 
to, and the cross-correlated spectral density for contacts 
n and m may be written in the formn3 



A d r V</>„ • V0 m • aT N (r). 



(18) 



Performing twice an integration by parts and making use 
of ( |l6| ) and a similar equation for / in T/v, equation (fill) 
is easily brought to a form 
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-4 / d 3 r<j) n m a / de(Vf) 



(19) 



As Tjv, G mn , and m are always positive, this proves 
that the current correlations are always negative, i. e. 
of fermionic type. As this proof holds for an arbitrary 
shape of the distribution function at the interfaces with 
electrodes, the correlations are negative for an arbitrary 
number of superconducting electrodes provided that the 
Josephson effect between them is suppressed. 

Consider now the noise at frequencies u> ~ eV. This is 
the range of quantum noise, which clearly cannot be de- 
scribed in terms of the Boltzmann equation. However 
it was recently shown that for disordered metals, the 
Langevin scheme may be extended to high frequenciesJay 
introducing a frequency-dependent noise temperatureEH 

T n {t,w) = ± Jde{f(T,e)[l-f(T,e + u)] 

+/(r,e + aO[l-/(r,e)]} (20) 

into Eq. (|l2]). If the frequency is smaller than 1/RC, 
where C is the capacity of the contact to all possible ex- 
ternal gates, the current fluctuation may still be consid- 
ered as constant along the contact length, and we arrive 
again at Eq. (13) with Tjv(x,w) substituted for Tjv(x). 
Note that the latter condition is well compatible with 
ui ^> D/L 2 because 1 jRC is larger by an extra factor con- 
taining the ratio of the contact diameter to the Thomas 
- Fermy screening length. 

Suppose first that eV < A. Integration in ( p0| ) and 
( |l3| ) then gives 



Si{w) 



1 

3~R 



[4U(w) + U(2eV + u) + U(2eV -w)], (21) 



with U(u>) = wcoth(w/2T), which reduces to 
Si(lo) = (l/3i?)[4w + max(4ey,2w)] 



(22) 



at zero temperature. The frequency dependence of Sj 
exhibits a kink at u) = 2eV, which would correspond 
to the Josephson frequency of the contact if the second 
electrode were also superconducting. The reason for the 
doubling of the kink frequency with respect to normal 
contacts is that the range of energies where the drop of 
the distrbution function takes place is two times larger 
now. 

In the case where T = and \eV\ > A, integration 
gives 

SiM = ¥ + i (|w + 2A| + |w - 2A| + \u - eV + A| 
K oil 

+\u - eV - A| + \lu + eV + A| + \lu + eV - A|) . (23) 

in accordance with recent numerical results,! this expres- 
sion exhibits peculiarities at oj = 2 A and lj = eV ± A, 
which are related with discontinuities in the energy de- 
pendence of the distribution function. The slope of S(u>) 
increases from 1/R at u) — ► to 2/R at li — > oo. 



In this work we have developped a semi-classical ap- 
proach for the noise in metallic diffusive NS-structures. 
This approach explains a number of key results in a phys- 
ically transparent manner and demonstrates that these 
results are immune to dephasing. 

This work was supported by the Swiss National Science 
Foundation. 
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